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Abstract 

We present a model description of the nucleon valence structure function 
applicable over the entire region of the Bjorken variable x, and above mod- 



ON ' erate values of Q 2 (~ 1 GeV 2 ). We stress the importance of describing the 

complete spectrum of intermediate states which are spectator to the deep- 
inelastic collision. At a scale of 1 GeV 2 the relevant degrees of freedom are 
constituent quarks and pions. The large- x region is then described in terms 
^ ' of scattering from constituent quarks in the nucleon, while the dressing of 



constituent quarks by pions plays an important role at intermediate x values. 
The correct small-x behavior, which is necessary for the proper normaliza- 
tion of the valence distributions, is guaranteed by modeling the asymptotic 

spectator mass spectrum according to Regge phenomenology. 
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I. INTRODUCTION 



Current deep inelastic scattering (DIS) experiments at HERA are probing the structure 
of the nucleon at ever smaller values of the Bjorken x variable (x = Q 2 /2Mu > 10~ 4 ), where 
M is the nucleon mass, and Q 2 and v the squared four-momentum and energy transfer to the 
nucleon. The large range of Q 2 available (4 < Q 2 < 3000 GeV 2 ) allows for quantitative tests 
of the nature of the evolution of structure functions at these low x [|I| . One reason why the 
small x region is interesting is that it is here that large distance phenomena become relevant, 
which offers the unique prospect of learning more about the nucleon's non-perturbative 
structure in DIS. 

Before we can fully appreciate the long range structure of the nucleon, there is of course 
the non-perturbative region at moderate x, but smaller Q 2 (Q 2 < 1 GeV 2 ), which has been 
the focus of many previous studies, but which still requires quantitative understanding. 
The significance of this kinematic region is that it reflects precisely the transition between 
the perturbative and non-perturbative domains of QCD, namely the interface between the 
parton picture of the nucleon, and the familiar valence quark models of hadrons which 
successfully describe much of the low-energy phenomenology. 

This connection has previously been investigated within phenomenological models of 
the nucleon, such as bag or non-relativistic quark models 0]. While these efforts have 
produced some encouraging results, some problems remain, however, in bridging the gap 
between leading-twist quark distributions calculated within such models, and the experi- 
mental structure function data. One problem encountered in most of these approaches has 
been their failure to correctly describe the x — ► behavior of valence quark distributions 
0. This is intrinsically related to the common approximation made in identifying the non- 
interacting quark system, which is formally on-mass-shell and remains spectator to the deep 
inelastic collision of the probe with the interacting quark, with a simple diquark, of mass 
iris ~ (2/3)M. A consequence of this assumption is that to ensure correct normalization, the 
calculated valence quark distribution must be evolved from extremely low resolution scales, 
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Q 2 ~ 0.1-0.3 GeV 2 0, where a s ~ 1-3. At these scales the use of perturbative QCD is 
rather questionable, and it's not even clear whether adding next-to- leading-order corrections 
is a sensible solution J7|. In this paper we will address the question of whether, and 
under what conditions, can the inclusion of higher mass spectator states generate sufficiently 
soft contributions to the valence quark densities so as to allow reliable descriptions of the 
data when evolved from a scale Q 2 ~ 1 GeV 2 , where perturbative QCD can be used with 
some more justification. 

Our physical framework is assumed to be as follows. At a scale Q 2 ~ 1 GeV 2 it will be 
natural to view the nucleon as composed of constituent quark (CQ) "quasi-particles" , which 
for our purposes will mean quarks having the quantum numbers of QCD (current) quarks, 
but which acquire a large average mass, rriq ~ M/3, due to their strong non-perturbative 
interactions. Since this scale is also typically that associated with chiral symmetry breaking, 
pseudoscalar Goldstone bosons (pions) will in addition play a role. We will therefore assume 
that the structure of the nucleon seen by a virtual photon with Q 2 ~ 1 GeV 2 is incorporated 
within the substructure of the CQs themselves and their pion dressing. 

In addition to knowing the internal structure of CQs, one also needs to describe the 
interaction between CQs in the nucleon, which we parametrize in terms of a relativistic 
quark-nucleon vertex function. We should state that our aim here is not to calculate the 
vertex function directly, and it will suffice to adopt a parametrization based on a physically 
motivated ansatz. Our main aim will be to model the spectral function which describes the 
non-interacting, spectator quark system, in terms of the degrees of freedom available to us 
at Q 2 ~ 1 GeV 2 , namely CQs and pions. Within our framework the task of describing the 
recoil spectator spectrum of the nucleon is reduced to detailing the remnant recoil system 
that remains after the interaction with the CQ. 

The outline of this paper is as follows. In Section [TT] we introduce our model for the spec- 
tral function of the spectator quark system. In Section |T| we discuss the CQ momentum 
distribution in the nucleon and the low mass spectator spectrum, which gives the leading 
contribution to the nucleon structure function at large x. The intermediate mass contribu- 
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tions, which are modeled by dressing the CQs by pions, are described in Section [TV|. Also 
discussed is the small- a; region of the quark distributions, for which we use a simple model 
based on Regge phenomenology to describe the large spectator mass continuum contribution 
to the spectral function. Finally, in Section [V] we discuss our numerical results and draw 
conclusions. 

II. SPECTRUM OF SPECTATOR STATES 

In this Section we outline the basic reasoning behind the need for a detailed description 
of the spectrum of states that are spectator to the interaction of the photon with a quark in 
the nucleon (see Fig.l). Let us begin the discussion with the familiar hadronic tensor W^ v 
which describes deep-inelastic scattering of a virtual photon (momentum q) from a nucleon 
target (momentum P): 

w»v(P,q) = ^J d ^ j q< (p\UOMo)\p) , (i) 

where J M is the electromagnetic current operator, and the states are normalized such that 
(P | P') = 2P (27r) 3 £ 4 (P — P'). For unpolarized scattering, which we consider throughout 
this paper, W^ u is usually decomposed into two independent structures: 

W,AP, 9) = [-9„ + 9 -f) F, + (p» - (p„ - q v ^j (2) 

where the structure functions Fig are expressed as functions of x and Q 2 {= — q 2 )- In the 
parton model Fig become independent of Q 2 (at leading order in a s ), and are related via: 
Fi(x) = F 2 (x)/2x = 1/2 J2q e q q{ x )i where q(x) gives the probability to find a quark q in 
the nucleon. 

In the usual treatment of inclusive DIS from a nucleon one formally sums over the 
complete set of intermediate states, which are labeled n in Fig.l, that are spectator to the 
deep-inelastic collision. In order to reliably calculate structure functions for all < x < 1 
one needs to know the spectator quark spectrum in its entirety, including the contributions 
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from the large spectator masses. In practice this is quite a challenge to any theoretical 
approach, and commonly one resorts to approximations which are, however, often valid only 
in a limited range of x 0. 

To elucidate the connection between the quark distribution q(x), calculated at some scale 
Q 2 = p 2 , and the spectator state n, we can use the so-called dispersion representation of the 
quark distribution |§||. Here q(x) is expressed in terms of the probability density, p, of the 
quark spectator states with invariant mass s = (P — p) 2 : 

2\ f fPmax( S ' X ) 2 / 2 2\ 

q{x;p)= ds dp p(s,p,x;n), (3) 

J J — oo 

where p is the four-momentum of the struck quark, and 

P 2 max (^x)=x( K M 2 - T ^j (4) 

is the kinematical maximum of the quark virtuality p 2 for given x and s. In terms of the 
amplitude ip n describing the absorption of a virtual photon by a quark with momentum p, 
this can be written as: 

p(s, p 2 , x; p 2 ) = J-^J-J2 ^n(P - P) i VV»(P - P) S(s - Ml), (5) 

where M 2 is the invariant mass of the intermediate spectator state n, and where formally 
the amplitude ip n (P — p) is defined via the matrix elements of the quark field operator at 
the origin ^(0): ip n (P — p) = (P — p, n \ip(0)\ P). An expression similar to (|3]) holds also 
for the antiquark distribution q(x). 

The dispersion representation in Eq.(|3p is valid as long as the integrals are convergent. 
Therefore it is generally assumed that the spectral density vanishes at large quark virtualities 
p 2 , and that the integration is dominated by contributions from the region of finite p 2 , 
\p 2 \ M 2 |||9|]. Then from Eq.(^j) one sees that the smaller the x, the broader the region of 
s from which q(x) receives contributions. For large x (x ^ 0.2), the region of finite s <i M 2 
is of major importance, and here one can expect models with low-mass spectators, such as 
diquarks, to be a good approximation. On the other hand, at small x the quark virtuality 
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p 2 is finite even for large s ~ M 2 /x. Therefore the quark distribution in Eq.([3D is sensitive 
to the high-energy part of the spectrum. 

The physical origin of the large-s region of the spectator quark spectrum is the high- 
energy scattering of the virtual qq component of the photon from the nucleon. This can 
be seen by recalling that the quark field operator acting on the target state can either 
annihilate a quark or create an antiquark in the target. In the former case the matrix 
element ip n oc (P — p, n\a(p)\P) describes the absorption of a virtual photon by quarks with 
momentum p, which results in a finite spectrum of intermediate states. On the other hand, 
the contribution from the antiquark part of the field operator, ip n oc (P — p, n\b\— p)\P), 
describes the "external" antiquarks from qq fluctuations of the virtual photon and their 
interaction with the nucleon. The possible antiquark momenta p are determined by the 
wave function of the photon and can be as large as the photon momentum q, leading to a 
large invariant mass of the antiquark-nucleon system, s = (P — p) 2 S> M 2 . 

The modeling of the spectral function p(s,p 2 , x; p 2 ) must of course reflect the resolution 
scale at which the target nucleon is probed. As mentioned in the Introduction, we will 
consider p at a fixed scale of order Q 2 = p 2 ~ 1 GeV 2 . At this model scale, constituent 
quarks are perhaps the most useful effective degrees of freedom. In a constituent quark 
picture, the resolution scale can be associated with some average size, Rq, of the CQ: 
p 2 ~ \jR 2 Q. Results for the color correlation function on a lattice ]T0|] indicate that color 



fields are confined to a distance of order 0.2 fm. This should also apply to each CQ, and one 
is led to the estimate p 2 ~ 1 GeV 2 . Naturally, at some other scale, CQs may not be the best 
variables with which to describe the spectral function. However, at p 2 ~ 1 GeV 2 , it is known 
from many other applications that CQs give a good description of the nucleon's dynamical 
properties, such as electromagnetic form factors. A natural question to ask therefore is 
whether CQs can be useful tools in describing DIS from the nucleon, Fig.2. 

Within this picture, it will be convenient to break up the spectrum of all possible spec- 
tator states into several mass regions, which we now summarize. 

• (.A): In a CQ picture, the simplest approximation is to identify the spectator system 
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which remains after one CQ is removed with a single diquark, with some relatively low mass, 
m s ~ 0(M), Fig. 3 (a). The low-s part (s < 1 GeV 2 ) of the spectral function can then be 
approximated by: 

p A (s,x,p 2 ) = 5 (s -m|) f Q / N (x } p 2 ) } (6) 

where the function /q/at(x,p 2 ), which is now independent of s, is identified with the distribu- 
tion of constituent quarks in the nucleon with a fraction of the nucleon light-cone momentum 
x and the invariant mass squared p 2 (see Eq.(^) below) (we omit explicit reference to the 
scale fi 2 ~ 1 GeV 2 from now on). In a relativistic treatment, p A is best described in terms of 
relativistic quark-diquark-nucleon (QDN) vertex functions. Simple effective models, with 
constituent quarks as basic ingredients, can be used to calculate the vertex function directly. 
Note that such calculations are only feasible if the complexity of the task is reduced to that 
of a few body problem. In other models, where the struck quark is treated as a current 
quark, the spectator state would have to be very complicated, making even an approximate 
model solution for the vertex function beyond our present technical ability. Therefore if the 
spectator state is to be treated as a diquark, then in order to solve for the relativistic vertex 
function, the struck quark must be a constituent quark. 

• (£>): To model the intermediate-mass part (1 GeV 2 ^5 s ^ sq) of the spectral function 
requires going beyond the simplest, single-spectator diquark approximation for the interme- 
diate state. One approach would be to consider higher-Fock state components containing 
a diquark plus some number of qq pairs as spectators. Within the bag model calculations 
of Ref. [PH| , for example, these component were modeled in terms of the insertion of an 



antiquark (or quark), produced from the incoming virtual photon, into the nucleon initial 
state, which then produced intermediate states with masses of order (5/4) M. 

At a scale Q 2 ~ 1 GeV 2 , as well as constituent quarks, pseudoscalar Goldstone bosons 
- pions — also enter as natural degrees of freedom in any effective model description. In a 
constituent quark picture, therefore, it will be more convenient to model these components 
in terms of the dressing of CQs by pions, as illustrated for example in Fig.3(b). In practice 



we restrict ourselves only the one-pion contributions. One could, in principle, continue 
the series in the pion number n to large n, hoping to describe the large-s region entirely in 
terms of pion ladder diagrams. However, taking this route would necessarily mean eventually 
including dressing due to the heavier meson spectators, as well as from diquarks. A more 
efficient approach will be to truncate the series at some mass s = Sq (which we take to be 
So ~ 2 GeV 2 ), and describe the large-mass contributions in terms of multi-quark and gluon 
configurations. 

• (C): To model the large (including asymptotic) mass tail of intermediate states (s > 
sq ~ 2 GeV 2 ), we describe the amplitude for scattering a qq pair from a CQ within Regge 
theory, which effectively takes into account an infinite set of ladders involving qq pairs in 
the intermediate state, Fig.3(c). The exchange of a Regge trajectory with the intercept 
«i?(~ 1/2) leads to imaginary parts of amplitudes rising with energy as E OLR 0. This 
energy dependence of the scattering amplitude then gives rise to valence quark distributions 
which rise at small x as x~ aR . 

Within the above framework, all three processes A-C will be described in terms of CQ 
parameters, with the basic common unit being the relativistic QDN vertex function. That 
is, the constituent quark-nucleon interaction is factored out from the hard scattering of the 
photon with the interacting quark. In the impulse approximation, nucleon DIS can then be 
viewed as a two-step process, in terms of photon-CQ, and CQ-nucleon subprocesses |12|-|14 



As we shall see, however, special care has to be given to the fact that the CQ is bound, and 
thus off-mass-shell. 



III. CONSTITUENT QUARK DISTRIBUTIONS IN THE NUCLEON 

At the relatively low momentum scales at which a constituent quark model may be 
relevant, higher twists (i.e. 1/Q 2 corrections to scaling quark distributions) are bound to 
play an important role. The final aim, however, is to compare the calculated structure 
functions with DIS data at large Q 2 (Q 2 > 5 — 10 GeV 2 ), where higher-twist effects should 



S 



be negligible. Therefore it will be sufficient to consider only the leading twist components 
of the structure function, represented by the "handbag" diagram in Fig.l. At leading twist, 
the fields that appear in the formal expressions for the current commutators in the forward 
virtual Compton scattering amplitude are those of point-like, elementary (i.e. current) 
quarks. In a constituent quark picture, however, the propagator of the quark before it 
is struck is modeled to be that of a quark with mass itlq ~ M/3 (while the propagating 
high-momentum quark after the interaction with the photon is a current quark). If one 
attempts to work with effective, composite degrees of freedom in DIS, a clear and well-defined 
connection must therefore exist between current and constituent quarks, in particular with 
respect to the nature of the photon-CQ vertex. 

An example of a simple dynamical model in which such a link exists is the Nambu 
& Jona-Lasinio (NJL) model [TIJ|n|. Here low- momentum gluons are integrated out and 
absorbed into a local, point-like interaction between quarks, characterized by an effective 
coupling constant. (Explicit perturbative gluonic degrees of freedom reappear when the 
calculated quark distributions are evolved to the higher Q 2 appropriate for the DIS region.) 
The mechanisms used in going from current to constituent quarks are dynamical quark mass 
generation, involving the scalar interaction of the quark with the Dirac vacuum, together 
with the spontaneous breaking of chiral symmetry. The connection between current and 



constituent quark masses is formally quantified via a "gap equation" ||15|| . Through inter- 
actions with quark condensates, (nearly) massless current quarks acquire large constituent 
masses, of order rug ~ 400 MeV. (For detailed reviews of the NJL approach see Refs. [|i~6|j.) 

Within this approach a CQ mass therefore appears in the quark propagator, while one 
still has a simple 7^ coupling at the photon-quark vertex. There is no form factor suppression 
at large Q 2 arising from the photon-CQ-current quark vertex. In a relativistic CQ picture, 
such as that embodied in the NJL model, the diagrams responsible for the CQ form factor 
are those in which the photon couples to the target CQ via qq loops. As Q 2 — > 00 these 
rapidly die out, leaving behind only the scaling part associated with the direct coupling 
of the photon to the CQ. In this way the only difference between scattering from current 
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and constituent quarks is therefore a renormalization of the quark mass. The dressing of 
the struck quark propagator is not included, since this would correspond to higher twist 
corrections. 



A. Quark Diquark Nucleon Vertex Functions 

Having made these remarks, we now turn to the calculation of constituent quark distri- 
butions inside the nucleon. Constituent quarks embedded in a nucleon have some charac- 
teristic momentum distribution, which is essentially given by the CQ wave function. Within 
the present approach it is more conveniently expressed in terms of the relativistic quark- 
diquark-nucleon vertex functions, Qqdn (proportional to the product of the wave function 
and the quark propagator), connecting the nucleon and off-shell quark to the non- interacting 
(spectator) intermediate state diquark. 

In the impulse approximation (i.e. no final state interactions between the struck quark 
and the residual diquark system) the CQ distribution (see Eq.®) can be written: 

Of.? 2 ) = / (v - £f ) * (f- - p 2 ) i {(p - p-r - mi) 

(7) 

where y is the fraction of light-cone momentum of the nucleon carried by the CQ, and the 
operator describes the soft CQ-nucleon interaction. The label (S) refers to the possible 
spin orientations of the spectator state (S — or 1), since any realistic model of the nucleon 
must incorporate pseudovector as well as scalar diquarks. In terms of the propagators and 
QDN vertex functions, TiS s > is defined as: 

n^ s=0 \P,p) = \{t-m Q )- X *% N (P,p) (P + M) SgUCP.p) ^-m Q y\ (8a) 
for a scalar (S = 0) spectator system, and: 

x(- 9 . + <^f^), < 8b > 
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OP . r, 



-Tr 



4rt s \p,p') 



for S = 1 spectators. Here mg is the CQ mass, and the conjugate vertex is defined as 

$QDN = 70 &QDN 70- 

The trace in Eq.(|7|) can be separated into terms proportional to the four- momenta P M 
and p^. 

^Tr[H^(P,p) 7/l ] = A (5 V) P, + /fV)P„, (9) 

where the functions and /2 s "* can De calculated from the QDN vertex. In terms of f^, 

(S) 

the constituent quark distribution function fq/ N is given by: 

Cd/^) = (/fV) + 1/ /f V)) ■ (io) 

Note that in general the functions depend on p 2 as well as P ■ p, however the latter is 
fixed by the 5- function in Eq.(^): P ■ p = (M 2 + p 2 — m|)/2. 

Different S = and S = 1 vertices, as well as a larger pseudovector diquark mass 
(7715=1 — ffi5 = o ~ 0.2 GeV), are needed to explicitly break spin-flavor SU(4) symmetry, 
which is reflected in a softer d quark distribution compared with the u quark distribution 
| |17| |. The complete Dirac structure for the vertex &qd N can ^ e written down in terms of a 
number of independent functions. In particular, one has for the S = vertex: 

^q > D n(P,p) = ^° ) I + < 0) (11) 
and similarly for the S" = 1 case: 

®q [ $n(P,p) = 7 a 7s + 0^ P a 75 + 0? P a 75 + • • • , (12) 

although different bases in which to expand &qdn can °f course be chosen fI5|| . 

In most previous calculations of structure functions involving QDN vertices [|T^-|23[ the 

( s) 

momentum dependence in Qqdn nas n °t been calculated directly, but rather has had to be 
parametrized, by appealing to DIS data to fix parameters. The approach adopted has been 
to choose a single form for both the S = and S = 1 vertices, and adjust the momentum 
dependence in each, to effectively compensate for the omission of the other structures. For 
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example, by choosing different shapes for the functions <fi( s=0 \ one can get similar results 
whether one uses a structure / fl9H2ll or one proportional to $ (see also Ref. ||23|| ). 

Although the phenomenological approach has been successful in allowing fair descriptions 
of the nucleon's valence quark distributions, it is naturally desirable to have stronger theo- 
retical justification for using any particular vertex structure. Recently attempts have been 

(S) 

made to calculate &qdn within effective models, notably the NJL model. Here we outline 



the general approaches used in these calculations — for more details see Refs. [I3|,f23-27 
The simplest ansatz has been to reduce the relativistic three-body calculation to a more 
tractable, effective two-body problem, by treating the nucleon as a composite of a quark 
and an elementary diquark |26| . Within this approximation one can solve exactly the Bethe- 



Salpeter equation for the quark-diquark system using a static approximation for the (local) 
quark exchange interaction. 

Beyond this simplest approximation, the Bethe-Salpeter equation has more recently been 
solved for the case where the exchanged quark was allowed to propagate between the quark 
and diquark fl8| , |2~8|| . Ishii et al. [^3] and Huang and Tjon [p5] have gone further by finding 



solutions to the covariant Faddeev equation without imposing the constraint that the diquark 



be elementary. Although only S = diquarks were modeled in Refs. p5| , |26[| , the authors of 
Refs. [0,0 have also incorporated pseudovector diquarks in their models for Qqdn- Note 
that to describe DIS one only needs the vertex function with an on-shell diquark — solutions 
to the Bethe-Salpeter or Faddeev equations more generally have both the quark and diquark 
off their mass shells. 

The results of the model calculations indicate that the S'-wave vertices give the 



most important contributions to the norm and to the nucleon's electric charge, which pro- 
vides some justification for the prescription applied in Refs. [[L9|-[23| of using only the leading 



structures. Furthermore, they show that the shapes of <p ( ' can be reasonably well approxi- 
mated by simple monopole or dipole functions: 



(m| - p 2 ) 



- (13) 
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with cut-off parameters As of the order of 1 — 2 GeV (depending on the precise value of ns). 
Such results are not surprising if one's model is to reproduce realistic values for the nucleon's 
static properties, such as the mass and r.m.s. radius. Note the presence of the zero in the 
numerator in Eq.(|i~3|) at p 2 = tuq, which has the role of canceling the quark propagators in 
fqw- m principle the quark propagators in Eq.([7|) could develop unphysical poles in the 
kinematically allowed region if the sum of quark and diquark masses falls below the nucleon 
mass, meaning that the nucleon could freely disintegrate into its constituents. As such the 
form factor in Eq.(^) reflects the underlying dynamics that are responsible for the bound 
state of the nucleon. Similar results have been reported for the quark-antiquark-meson 
vertex function in Ref. p9 . 



Note that As is not related to the intrinsic NJL cut-off, Ajvjl- By treating scalar or 
pseudovector diquarks as quasi-particles, with effective boson propagators oc ((P — p) 2 — 



" 1 , one renders all loop integrals appearing in the calculation of the Bethe-Salpeter 



amplitudes finite, without any need to introduce ultra-violet cut-offs as regulators. This can 
be easily seen by counting inverse powers of momentum running through the quark-diquark 
legs. The cut-off As is therefore only an indicator of the "hardness" of the QDN vertex, or 
the "size" of the quark-diquark configuration in the nucleon. 

For the scalar vertex ([TT|), the S'-wave structure <£> { q ] dn = I <f><® |U||0] in Eq.(0) gives 
for the functions f^: 



fl°\p 2 ) 
/ 2 ( °V) 



2 (m 2 Q - p 2 ) ' 



{m Q M + P-p) 



(14a) 
(14b) 



(m 2 Q - p 

Note that an equivalent result is obtained from the structure P, or (/+ P/M)/2, which 
would be the leading one in the non-relativistic limit With the pseudovector vertex 



structure: $ 



(i) 

QDN 



757«0^ one has: 



(1 V) 2 (p 2 -M 2 -2m 2 ) 



2 (p 2 - m 2 Q ) 



nit 



(15a) 
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/ 2 (1 V) = ~ - ' \ v " v r v - m 



(1) (P 2 ) 2 /(m 2 -M 2 ) 



2 y\#2^ 



2 (p 2 - 



2 \2 



mt 



+ M 2 + 6Mm Q - m 2 Q - 2m 2 s + 2p 2 j . (15b) 

From Eqs.([|), (|7|) - (|1|) we can now calculate the quark distribution associated with 
process A, which, because of the low mass present in the spectator state, should be a good 
approximation at large values of x. 

B. Nucleon Structure Functions at Large x 

At large x one can think of the photon scattering quasi-elastically from the constituent 
quark, in the sense that the internal structure of the CQ "quasi-particle" is not resolved. 
In this case the CQ structure function is a 5-function, and the resulting nucleon quark 
distribution given entirely by the momentum distributions fQ^iViP 2 )- Integrating fqj N 
over p 2 gives then the valence quark distribution in the proton: 

Qa ( X )= — / d P fQ/NW'P ) 6 ( 1 ~ X Q) = / d P fQN\. X >V)> ( 16 ) 

Jx y J— oo J— CO 

where x = Q 2 /2P ■ q, xq = x/y is the Bjorken variable of the CQ and the maximum value 
of the quark virtuality is given by Eq.(|j). 

The parameters defining the constituent quark distribution fg) N (and hence q^) are as 
follows. Consistent with typical masses in effective constituent quark models, we take for the 
CQ mass uiq w 450 MeV. For the scalar (pseudovector) diquark mass we use mom ~ 1.0(1.2) 
GeV. These diquark masses are somewhat larger than in some model calculations of the 
nucleon's static properties, which reflects the approximation that the low-mass spectator 
system can be described as a single diquark, Eq.(^). In fact, one would expect in effective 



CQ-diquark models to see a spectrum of diquark bound states with a range of masses Jib 
By taking into account only the lightest diquark state one necessarily underestimates the 
average mass of the diquark which would effectively describe the low-s part of the spectrum. 
For comparison, we also examine the effect on the quark distributions of using smaller 
masses, m (i) ~ 0.6(0.8) GeV. 
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For the QDN form factor (vertex function) specifying the distributions fgo^ we use 
the form in Eq.(|l3|). The exponential parameters ns in Eq . (|13D are fixed by the behavior of 
the quark distribution at large x, since 

From Eqs. ([l7|) and the expressions for fq) N one finds that with no = 2.0 and n\ = 3.5 the 
quark distributions q ( -°' 1 \x) — > (1 — x) 3,4 as x — > 1 [30]. The values of the cut-off parameters 



As in Eq.(|l3|) are taken to be: A = 1.0 GeV for the scalar vertex, and Ai = 1.2 GeV for 
the 757„ vertex. 

Using these parameters, we show in Fig. 4(a) the valence sum (solid curve), x(uy + dy) — > 
3x (qjp + /2, evolved from u 2 = 1 GeV 2 to Q 2 = 10 GeV 2 . One sees that within the 
impulse approximation the valence quark distribution is almost a factor 2 larger than the 
data [^OJ (shaded) at Q 2 = 10 GeV 2 . The position of the peak value of x{uy + dy) depends 
quite sensitively on the mass of the intermediate state diquark system. For the smaller 
diquark masses (mom = 0.6(0.8) GeV), the result (dashed) is nearly a factor 3 larger, and 
peaks at a considerably larger value of x compared with the data. The valence quark ratio 
in the impulse approximation, dy/uy — > 2q^ I (Sq^ H-^), shown in Fig.4(b), exhibits 
far less dependence on the diquark masses, and gives quite reasonable agreement with the 
data at large x, particularly for the larger mass spectators. 

The fact that the impulse approximation is only good at large x is clearly reflected in 
the absence of the correct (singular) behavior of q^\x) as x — > that would be expected 
from Regge phenomenology, namely ~ x~ aR (an ~ 1/2) (see Section IV.B). Because only 
small (finite) spectator masses are included in q^\x), the calculated distributions are finite 
at x = 0, whereas the experimental distributions diverge. Therefore a large fraction of the 
normalization has to come from the medium and large-x regions to compensate — hence 
the large peak in x(uy + dy). This fact demonstrates that the distributions q A (x) alone, 
calculated in a CQ model, are incapable of describing the valence structure function over 
the entire range of x, and points to the necessity of including larger spectator masses, 
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corresponding to processes B and C outlined in Section II. In the next Section we describe 
the contributions to the structure functions from the larger spectator mass components of 
the quark spectral function. 



IV. STRUCTURE OF CONSTITUENT QUARKS 

As clearly seen in the previous Section, one should not expect valence quark distributions 
calculated from the handbag diagram alone to exhibit the correct behavior at small x. Since 
this is intrinsically related to the presence of low-mass spectators in the intermediate state, 
to rectify this deficiency one must consider in addition contributions to the spectral function 
p from states with masses larger than irts ~ M. 

A systematic formulation within the CQ picture developed in Section III is one in which 
the interactions among quarks are completely factored out of the basic QDN vertex. Inter- 
actions between quarks then give rise to effective constituent quark dressing — for process 
B by a cloud of pions in the intermediate- s region, while for process C by a soft "Reggeon 
cloud" at large s — which would result in CQs attaining substructure, Fig.2. The contri- 
bution from these processes to the quark distribution (direct photon diagram, Fig. 2 (a)) can 
then be written: 

d A p 



<&%) = J j^2ir5((P-pf-ml)Tr[Q(p, q )H( s \P,p)}, (18) 
where the operator 



Q(P, 0) = ~i J S ((k + qf) [4{Jk- niQ)- 1 T(p, k)(fe- mg)" 1 ] , (19) 

describes DIS from a dressed (generally off-shell) constituent quark. The trace in Eq.fll8|) 
is taken over the upper and lower quark indices in Fig.2 (a). The crossed photon diagram 
in Fig.2 (b) gives rise to an antiquark distribution q(x), which is given by an expression 
analogous to Eq. (|19|) with the substitution q — > —q. 

The operator T(p,k) in Eq.(|I^) describes the truncated quark-quark interaction. In 
what follows we assume that T is an analytic function of the invariant mass squared of 
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the quark-quark system, sq = (p — k) 2 , of uq = (p + k) 2 , and of the quark virtualities 
k 2 and p 2 . For real sq and uq, T has a right-hand cut in the variable sq, a left-hand cut 
in uq, and singularities for p 2 > and A; 2 > 0. To make use of these properties in the 
loop integral in Eq. (|i9|) , it is convenient to parametrize the loop momentum k in terms of 
external momenta p and q: k = ap + f3q' + kx, with q' = xqP + q, Xq = —q 2 /2p ■ q being the 
Bjorken variable of the CQ, and kr is a two-dimensional vector perpendicular to both p and 
q. Then the 5-function in ( |19|) fixes a = xq. After integrating with respect to (3, one finds 
that the operator Q vanishes outside the interval \xq\ < 1. For < xq < 1, Q is given by 
a dispersion integral in sq along the right-hand (R) cut, 

11 r d 2 k T ds Q \m R 4(}k + mQ)T(p,k)(}i + mQ) 
q) -TQ^2p—qJ T^x- Q (F-m 2 ,) 2 ' (2 ° } 

and hence contributes to the quark distribution. For —1 < xq < 0, the operator Q is given 
by the dispersion integral along the (left-hand, L) cut in the w-channel, 

Q , v _J 1_ r d 2 k T du Q Im L + m Q )T(p, fc)(# + m Q ) 

^LlP,gj iQ^ 2 p-qJ l+x Q (k 2 -m 2 Q ) 2 ' 1 J 

determines in fact the contribution from the crossed photon diagram in the physical 
region xq > 0, and corresponds to the antiquark distribution in the CQ. To obtain the 
antiquark operator Q from Eq.(^) one must substitute q — ► — g (and therefore Xq — > — xq) 
and introduce an overall minus sign, Q(p, q) = —Ql(p, —<l) (see e.g. 0). 

Before discussing the details of the quark-quark interaction within our model, we can 
make one more model-independent observation about factorization of the total 7*iV inter- 
action. Let us firstly expand the operator Q in terms of Dirac basis tensors: 

G ^ = I ^ + >^ + '^i- (22) 

where the coefficients Qi (i = 0, 1, 2) scalar dimensionless functions of xq and p 2 . An 
expression similar to ( |22|) holds for the antiquark operator with the scalar coefficients <2o,i,2- 
Projecting Q onto positive energy space and averaging over CQ spins one can define the 
quark and antiquark distributions in a CQ as: 

17 



^Tr [(y + m Q ) Q{p, q)} = Q + ^ T Q 1 + Q 2 , 
2 " rrin 



-Tr 

2 



(t + m Q )Q(p,q) = Q Q + J^Q 1 + Q 2 . 



(23a) 
(23b) 



For a quark (antiquark) with no internal structure one can easily show that in the Bjorken 
limit only a single function is non-zero: Q 2 ^ (<2 2 7^ 0), while <2 ,i = (Qo,i = 0). On 
the other hand, for a constituent quark with internal structure all three coefficients can be 
non-zero and therefore contribute to the quark (antiquark) distribution. 

Finally, substituting (^) into Eq.([L8|) and using Eqs.(|9]), (|T0D and ( |23a| ), the contribution 
to the nucleon's quark distribution from dressed CQs is: 



x y 



dp 2 



fQ/N(y>p' 



X 



JJ 



,P 



(24) 



where fQ S J N {y,p 2 ) is given by fllPf) while the functions A/q^ are 

&ft> s \y,P 2 ) = ■^T Tc {rt s) (P,p)} - f ( Q %(y,P 2 ), 

a AS)/ 2n yP-p-p 2 f (S), 
S7T 77tr> 



(25a) 
(25b) 



A similar expression can be written for the antiquark distribution q(x) with the substitution 
Qi -> Qi and q Q -> g Q . 

The first term in Eq.fl2~4|) is a generalization of Eq . (|16D for an 'undressed' CQ to the 
case of a CQ having non-trivial internal structure, and is written as a convolution of the 



CQ distribution fh) N and the spin averaged CQ structure function qg, Eq.( [23a|) . One sees, 



however, that the non-trivial Dirac structure of the 7*Q interaction leads to the convolution 
breaking terms in Eq.( 



which thus render the convolution hypothesis p3| , ^4| , |36| invalid. 
Note also that a simple convolution formulation of the structure functions is recovered in 

the p 2 — > uiq limit, since the terms A/0,1 come with an extra factor of (p 2 — rug) compared 

ts) 

with the fq/ N term. 

As discussed in Section [D], it is convenient to split the spectrum of invariant masses of 
the quark-quark system into two regions: small Sq (sq < M 2 ) and large Sq (sq > M 2 ). In 



what follows we consider a model where the interaction T is described by meson exchanges 
for the low-mass region, while to evaluate the effect from the large-mass region we adopt a 
simple model based on Regge phenomenology. 



A. Pion-Dressed Constituent Quarks (Process £>) 

In evaluating the effect of the dressing of CQs by meson "clouds" , we will restrict our- 
selves to pions. Contributions from CQs dressed by higher mass mesons (e.g. p, u, etc.) 
will be largely suppressed by the larger meson masses |$4|, and only ever become noticeable 



when unrealistically hard form factors are employed. In addition, predictions for the p and 
heavier mesons are less reliable in constituent quark models due to threshold effects, if one 
uses CQ masses of the order ~ 400 MeV. 

For the interaction of CQs with pions we will utilize the simple effective chiral Lagrangian, 
valid at scales below ~ 1 GeV, similar to that considered in Refs. p5|j36[1 , in which the leading 
interaction is given by: 

£qtt = ^5 da ^- ( 26 ) 

Here ipQ and <j) n are the CQ and pion fields, g& ~ 0.75 is the axial coupling constant of 
the constituent quark, and = 93 MeV is the pion decay constant. A typical process 
which arises from this interaction Lagrangian is depicted in Fig. 5. Note that the pion fields 
appearing in Cq-k represent point-like pions. In principle Cq-x gives rise therefore also to 
terms involving the direct coupling of photons to pions. At high energies, however, the 
pion's substructure would be expected to play a role, and these contributions, through the 
presence of the pion form factor, would be suppressed. Furthermore, the direct 7*7r couplings 
would also violate the Callan-Gross relation, and as such they are not considered. 

In the present work we consider only order 0(1/ f n ) effects from Cq^. The relevant 
process to be considered is therefore the one-pion exchange process in Fig. 6 (together with 
the crossed photon diagram). Note that scattering from the pion cloud itself, as would be 
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needed to generate antiquark distributions, for example, would enter only at order 0(1/ f^) 2 
from the Lagrangian |J7|. The pion cloud contributions were considered in Ref. ||36|| , however, 
the expansion of Cqt, was performed only up to 0(1/ f^). 

The quark-quark operator T for the one pion exchange interaction can then be obtained 
from fl2"S| ) by a Fierz transformation: 

lm R T* = 2n5 ((p - k) 2 - m 2 ) T (r* T) , (27) 

where T n = (g^/ \plj. K )(p' — ^)15<Pq-k represents the QQtt vertex, with (pQ n being the QQn 
vertex function, and T is the Fierz transformation operator. Substituting fl27|) into (^0|), we 
find the one pion exchange contribution to the coefficient functions Q,{. 

Jl " 8^ V2U ) loo (k 2 - m 2 Q ) 2 

x {(1 - x Q ) (m 2 Q -k 2 + x Q (p 2 - m 2 Q )) - x Q m 2 } , (28b) 

J_ / 2 2 0^ 1 

^ 2 8tt» ^ ^ ; loo (P - m 2 Q 

x { (p 2 (2x Q -l)-k 2 ) (m 2 Q - k 2 + x Q (p 2 - m 2 Q )) + x Q m\ (p 2 - m 2 Q ) } . (28c) 

The kinematical maximum of the quark virtuality, k^ ax (sQ, xq,p 2 ), is given by Eq.(f|), with 
sq = m 2 and the nucleon mass squared M 2 replaced by p 2 . 

Equations (|28|) determine the pion contribution to the quark distribution in the CQ via 
Eq. (p4[) . Similar expressions arise for the pion contribution to the antiquark distribution, 
~Q* , and are related to (|28"D by crossing symmetry, viz. interchanging xq — > —Xq and 
introducing an overall minus sign. In the present work we consider the pion contribution to 
the valence part of the nucleon quark distribution, so that only terms even in xq survive. 
In this case the scalar term Qq cancels out and the convolution breaking term in Eq.(p4|) 
arises only from the Dirac structure p 1 in Eq. (p2|) . 

In Eqs.(|2"E|) the integral over k 2 is regularized by a sharp cut-off, Aq^, whereby the QQn 
vertex function is parametrized by a theta-function: 
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Mk\p*) = Q (k 2 + A 2 Qt ) 6 (p 2 + A 2 QlT ) , (29) 

where we have also parametrized the p 2 dependence of this vertex by a similar cut-off (for 
other forms of the cut-off see Refs. p6|j38|j39[| ) . The parameter Aq^ represents a typical scale 



beyond which the effective interaction in Eq.flzq) is no longer a valid approximation to QCD, 
and is usually taken to be of the order Aq w ~ 1 GeV P5fl . 

In Fig. 7 we plot the one-pion exchange contribution to the sum of the valence 
u v (- 3/2 ( q f - + 5/2 (g£> - $>)) and d v (- 3 - + 2 ( q $> - $>)) 
distributions, compared with the impulse approximation contribution, for three values of 
the ultra-violet cut-off Ag n , namely 0.8, 1.0 and 1.2 GeV. The results indicate that for 
reasonable values of Aq w (< 1 GeV), the one-pion exchange contribution is certainly not 
negligible for x ^ 0.2, and indeed provides around 15-20% of the valence quark normalization 
compared to the impulse approximation distributions associated with process A. (Similar 
numbers are obtained in non-covariant, infinite momentum frame calculations, where one 
makes use of a transverse momentum cut-off, (j)Q w = O (k^ — Agjj — see Appendix A.) 
However, it is also apparent that DIS from dressed CQ is not able to provide sufficiently 
soft contributions to the nucleon quark distribution that would simulate Regge behavior. 
Agreement with large Q 2 DIS data would still require evolution from scales similar to that 
needed for process A alone, fi 2 ~ 0.2 GeV 2 p3 |. Adding higher order contributions would 



make the total distribution softer, however, one would need to include quite high orders 
to come close to reproducing the required degree of singularity as x — > 0. The complexity 
involved in performing even the 2-loop calculation |37j based on the effective interaction 
Cqk is quite daunting! In the next Section we shall examine a more efficient method of 
incorporating the larger mass continuum of spectator states (above s ~ 2 GeV 2 ), based on 
Regge phenomenology. 
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B. Asymptotic Spectator Masses (Process C) 

It should be clear by now that to obtain quark distributions which are sufficiently soft 
at small x to be compatible with DIS data (when evolved from our model scale of /i 2 ~ 1 
GeV 2 ) one must include contributions to the spectator spectral function from the large-mass 
continuum. As mentioned earlier (see also Ref. ||), the physical mechanism to which one 
can attribute the large mass states is the high-energy scattering of the qq components of 
the photon from the nucleon. By examining the structure of the energy denominators of 
the qq + N system in the target rest frame one can immediately see that this time-ordering 
dominates the contribution from the direct scattering of the photon from a quark in the 
target. The virtual qq fluctuation excites the target into states with a spectrum of masses 
s > 0(M 2 ). 

In modeling the large-s contribution to the spectral function it will be convenient to 
make use of well-known ideas from the Regge poles approach to high-energy scattering ||. 
Formally, the sum over the spectator masses s can be shown to be effectively described by the 
exchange of Regge trajectories in the complex angular momentum plane. In a constituent 
quark model, where the interaction of the photon with the CQ is factored from the total 
amplitude, the qq pair will in fact excite the CQ in the nucleon, which necessitates integrating 
over states that are spectator to the collision. The problem then reduces to describing 
the asymptotic behavior of the quark-quark interaction T(p, k) in terms of exchange of 
Regge poles, Fig. 3(c) (for an alternative description of high-energy quark-quark scattering 
see Ref. @.) 

The leading contribution to the asymptotic scattering amplitude comes from the Pomeron 
exchange, which, however, cancels out in the valence quark distribution. The relevant Regge 
trajectories therefore correspond to the exchange of vector mesons, and we will assume that 
the corresponding Reggeons couple to quarks like vector mesons. In this case the quark- 
quark scattering amplitude at high energies can be written: 
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T{p, k) = 1{q)tx ^ , (30) 

where T is a Lorentz-invariant function, whose interpretation becomes clear if we project T 
onto positive energy space and average over quark spins: 

T(s Q , u Q] k\ p 2 ) = ^Tr 9 Tr [(- # + m Q ) {p 1 + m Q ) T(p, k)} . (31) 

Here the trace is taken over the upper and lower quark indices in Fig. 2, which we distinguish 
by indices (g) and (Q), respectively. Therefore the function T(sq,uq] k 2 ,p 2 ) can be consid- 
ered as an analytical continuation of the spin-averaged quark-quark scattering amplitude 
into the off-mass- shell region. In general T describes both the qQ and qQ scattering chan- 
nels, which are connected by crossing symmetry. Namely, in the s-channel T(sq,uq; k 2 ,p 2 ) 
coincides with the qQ scattering amplitude, T 5 q(sq; k 2 ,p 2 ), which contributes to the quark 
distribution, while in the w-channel it is identified with the qQ amplitude, T 9 q(mq; k 2 ,p 2 ), 
from which one obtains the antiquark distribution. 

For the interaction (|3~0"D , the coefficients Q§ x 2 in the high-energy (Regge) region become: 



Qo = 0, (32a) 

ef ~ [ds Q f k - dk 2 frg«(«*y) f 1 1 (32b) 



(2n) 3 J y-oo (k 2 -m 2 Q ) \s Q -p 2 -k 2 

y*_ 1 f, [ k ™* , l2 lmTg Q {s Q ;k 2 ,p 2 ) \ m 2 Q - k 2 - 2x 2 Q p 2 



= 77^ / ds Q / dk 2 gvv V1 ' -x Q + 2 X > ( 32c ) 

(2tt)3 7 y 7_oo (k 2 -m 2 Q ) I sq-p 2 -A; 2 J 

where fc^ aa , = fc^ aa .(sQ, xq,p 2 ) is given by Eq.(^), with the nucleon mass squared M 2 replaced 
by p 2 . Together with Eq.([24|), these give the contribution to the nucleon quark distribution 
from the qq scattering mechanism. Expressions for the antiquark distributions, Q 1 2 , can 
be obtained by applying the crossing symmetry rules. The final result is similar to Eqs.(j32j), 
with the qQ scattering amplitude T 5 q replaced by the qQ scattering amplitude T q Q. 

To finally calculate the Regge-exchange contribution to the quark (qt) and antiquark 
(q~c) distributions from Eqs. (|24| , |3"2"|) requires modeling the qQ and qQ scattering amplitudes. 
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As discussed in Ref. 0, at high energy and low x one can approximate TqQ and T q Q by 
the constituent quark-target amplitudes obtained from the simple additive quark model of 
high-energy hadron scattering. In this case, the nucleon-nucleon scattering amplitude can 
be expressed in terms of the quark-quark amplitude and the CQ distribution functions /q/jv : 



= / [dp] [dp] f Q /N [y,p ) fq/N [y,p J (33) 

sn J \ / \ / S q 

where is now the total hadronic center of mass energy squared, p and p' are the four- 
momenta of CQs in nucleons iV and N', with y and y' the light-cone momentum fractions 
carried by quarks Q and Q', and sq = (p+p') 2 is the invariant mass squared of the interacting 
QQ' system. The brackets in Eq. (|33"D denote: J[dp] = §q dy/y j 1 ^^™ 3 ^ dp 2 . At high 
energy, where sq ~ y y' s^, one can obtain approximate solutions to Eq. (|33|) by evaluating 
the quark-quark amplitude at averaged values of y and p 2 . With the QDN vertex model 
parameters in Section III, one finds (y) ~ 1/3 and (p 2 ) ~ —(0.3 — 0.4) GeV 2 . 

The SQ-dependence of the amplitudes corresponding to the valence quark distribution, 
namely the difference ATqq = T^q — T q q, is then fixed directly by the energy dependence 
of the pp and pp amplitude difference. In the Regge approach this difference is determined 
by the exchange of the spin-1 omeson trajectory: 

Im (T Pp (s N ) - T pp (s N )) = R w , (34) 

where is the residue of the uj Reggeon trajectory, and the intercept. The best fit of 



Ref. |42| gives R w ~ 42 mb • GeV 1 a " and a u = 0.548 . For the p 2 and k 2 dependence of 



ATqq we follow Ref. || in adopting the following factorized form: 

AT QQ (s Q ;k 2 ,p 2 ) = g R (k 2 ) g R (p 2 ) AT QQ (s Q ; 0, 0), (35) 

where the form factor functions g R describe the fall-off with four-momentum of the total 
quark-quark amplitude: 

^ = (1-^. (36) 
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For the exponent ur we choose the value ur ~ 4 to correctly reproduce the tail of the 
resulting distribution function at large x, while the cut-off Ar remains a free parameter, to 
be determined by fitting to the DIS data. The best fit to the data gives Ar w 0.25 GeV. 

Finally, with these ingredients the contributions to the uy and dy valence quark distri- 
butions from process C are both given by 3/2 ((<7c°^ ~ Se ) + ylc 1 ~ 9e ))> where qf^ are 
obtained from Eqs. fl24]) , (|25|) and and (p2]). 



V. RESULTS AND CONCLUSIONS 

With the inclusion of all three processes *4, B and C, describing the low, intermediate and 
high-mass spectrum of spectator states, respectively, the valence u and d quark distribution 
are given by: 

u v {x) = Z u {u A (x) + (u B (x) -u B (x)) + (u c (x) - u c (x))} , (37a) 
d v (x) = Z d \d A (x) + (d B (x) - d B (x)^ + (dc(x) - dc(xf\\ . (37b) 



The normalization constants Z U ( i are determined by charge conservation, 

2 



2 + (u — u)b + (u — u) 



c 



(38a) 



Zd ~ i + {d _d) B + (d- d) c ' (38b) 

where the brackets (• • ■) denote the first moment, and where the impulse approximation 
contributions, d A and u A are normalized to 1 and 2, respectively. 

The relative normalizations of the individual contributions from A — C are completely 
fixed in terms of the parameters described in Sections III and IV. The contributions to the 
normalization of the valence uy + dy distribution from one pion exchange is ~ 16%, while 
the Regge exchange contribution is ~ 28%. This means that the impulse approximation 
contribution (process A) is renormalized by Z u = 62% and Zd = 43% for the u and d 
distributions, respectively, and contributes about 56% to the total normalization. The con- 
tributions to the second moments of uy + dy at the input scale are ~ 41%, 6% and 0.06%, 
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so that some 53% of the nucleoli's momentum at /i 2 = 1 GeV 2 is left to be carried by sea 
quarks and gluons. 

The results for x(uv + dv), evolved from /i 2 = 1 GeV 2 , are plotted in Fig.8, in comparison 
with the data (shaded) at Q 2 = 10 GeV 2 . Also shown for comparison are the corresponding 
results for processes A and £>, and A alone. In total, the results indicate that one can indeed 
obtain a good description of data above /z 2 = 1 GeV 2 in terms of CQ parameters. This scale, 
being at the boundary of the scales which are generally accepted as reliable for perturbative 
evolution, is considerably larger than those used in many previous attempts to calculate 
leading twist quark distributions (which typically have needed /i 2 ~ 0.1 — 0.3 GeV 2 ). This 
is one of the main results of the present work. 

As possible extensions of this work, one could next attempt to incorporate sea quark 
distributions into the formalism at small x, which would require a more detailed analysis of 
processes B and C. In the former, the sea would simply be associated with the antiquark 
content of the pion cloud. However, to adequately describe the low-x region one would 
need to extend the Reggeon model for process C to include in addition the exchange of 
Pomeron trajectories between quarks. Another possible test of the model would be to 
consider polarization observables, such as the nucleon's spin-dependent g\ structure function. 

Finally, we should state that while our formalism is quite general, it is nevertheless a 
formidable challenge to embed within any specific model the full spectator mass spectrum, so 
as to describe consistently the nucleon's structure function over the whole range of Bjorken- 
x. Our attempt here to obtain such a unified description suggests, however, that constituent 
quarks may still be the most relevant degrees of freedom appropriate for this task, and for 
the more general problem of understanding the nucleon's non-perturbative structure in DIS. 
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APPENDIX A: CONSTITUENT QUARKS IN THE INFINITE MOMENTUM 

FRAME 



In this Appendix we discuss an alternative, non-covariant, formulation of DIS from con- 
stituent quarks to that in Section IV. In a covariant formalism contributions to the quark 
distribution functions from DIS off dressed CQs can be expressed as a two-dimensional con- 
volution (in y and p 2 ) of the CQ distribution fq/ N and the off-shell CQ structure function 
q<2, plus a non-convolution term whose contribution is proportional to (p 2 — itlq). In the on- 
mass-shell limit this term would vanish, leaving only the convolution term. Furthermore, the 
remaining term would reduce to a one- dimensional convolution in the momentum fraction 
y only. From a practical point of view, this would considerably reduce the number of inte- 
grations needed, and simplify calculations involving higher-order, multi-loop, pion-dressing 
contributions. 

A framework in which one-dimensional convolution equations can be recovered is non- 
covariant time-ordered perturbation theory in the infinite momentum frame (IMF). Because 
here all particles are on-mass-shell, the non-convolution off-shell corrections in Eqs.(|24]) are 
identically zero, and factorization of subprocesses is automatic. Furthermore, in the IMF the 
so-called Z-diagrams, which involve antiparticles, or particles moving backwards in time, are 
suppressed by powers of 1/Pl P2| , where P L is the longitudinal momentum of the target. 
Hence only forward moving particles are ever considered. It is only within this approach 
that one can consider the function fqw (integrated over transverse momentum) as a genuine 
probability distribution function — in any other frame one must also incorporate antiquark 
components into the virtual constituent quark (the IMF is also the most appropriate context 
in which to view the q — > qn splitting functions of Eichten et al. and the nucleon — > 

nucleon + tc distributions of Refs. p3| , p4| .) Therefore the simple one- dimensional convolution 
formulation is necessarily frame-dependent. 

By way of illustration, let us examine the one-pion exchange contribution to the quark 
distribution function, as in Fig. 6, in the IMF. The momentum of a nucleon moving with 
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longitudinal momentum Pl in the IMF can be parametrized as f32U33[ 



M 2 

Pl + 2PZ ;0t ' Pl 



Similarly for the CQ and spectator diquark momenta, we have: 

p,= [\y\P L + „?_.,„ ; pt, yP L 



(Ala) 



P. 



2\v\Pl 



2|l-y|P£ 

and for the struck quark and spectator pion: 



(Alb) 
(Ale) 



ki, 



1 \ ID , ™Q + ( k T + fPT 

fI-PlH 

V 



k T + -p T , xP L 



k' = \y-x\P L + 



2\x\P L 

ml + (-k T + (l - |) p T 

2|y-x|P L 

-k T + h --Jpr, (y-x)P L 



(Aid) 



(Ale) 



respectively. Evaluating the trace as in Eq.(|19D in terms of these momenta, one finds the 
contribution to the valence quark distribution from a constituent quark dressed by a pion: 



q ( i\x) 



ri dy roc . . poo 

— dp 2 T f ( Q L(y,p T ) / dk 2 T q Qn/Q (x Q ,k T ), 
Jx y Jo jo 



(A2) 



where in the IMF the functions fn N are given by 



Q/N 



fQ/ N (y,PT) 



for S = spectators, and 



f ( Q/ N (yiPT) 



<f>W{y, PT ) 



16tt 2 y 2 (M 2 - M 2 ) 2 



| p 2 + ( mQ+yM ) 2 } 



(A3) 



16tt 2 y (M 2 - M 2 QD ) 2 
for 5 = 1 intermediate states. The variable 



6Mm Q + 2P-p + 



AP ■ p' p ■ p' 



mi 



(A4) 



2 rriq+PT m 2 s +p 2 T 



MA D = ( P +p'y 



y 



i-y 



(A5) 
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is the squared mass of the virtual quark-spectator diquark system. Note that the IMF 
expressions for fqJ N {y,PT) can be related to the covariant expressions in Section III (apart 
from the QDN vertex function) if one observes that p^ in Eq.( [Alb| ) satisfies: 

2 

where p 2 max is given by Eq.(|). 
The function: 

r 7 \ 1 ( 9Am Q \ 2 \<p QlT {x Q ,kT)\ 2 r ;2 . 2 /-, n2\ / Av n 

represents the (unintegrated) structure function of a CQ dressed by a pion (see Eqs.( [23a| ), 
(f2~8|)), where the squared mass of the Qir state is given by: 



M^(* + *f=M + lM (A8) 



One can easily demonstrate that a pseudoscalar Qn interaction leads to the same result as 
in Eq. ([A7|) , by virtue of the IMF on- mass-shell condition for the CQs (i.e. = mq in 

Eq.(pHD), and the Goldberger-Treiman relation fT6|| . 

Although the kinematical (or trace) factors are similar in the covariant and IMF calcu- 
lations, the connection between the vertex functions (or form factors) in the two approaches 
is not a straightforward one. Covariantly, the vertex function such as in Eq.(^) can only 
depend on the virtuality of the exchanged off-mass-shell particle, p 2 . In the time-ordered 
IMF approach the vertex function cannot depend on the p 2 , but rather must be a func- 
tion of the variable Mq D . This constraint stems from the requirement that the momentum 
distribution functions respect probability and momentum conservation, hence must display 
explicit symmetry under the interchange y <-> 1 — y (Eq.( [Bl|) in Appendix B below), which 
cannot be the case in a covariant formulation. 



Since the pt and kx dependence in Eq . (|A2|) is factorized, one can write the total quark 
distribution as a one- dimensional convolution of the CQ structure function, qq, with the 
Pr-integrated CQ distribution fg/ N (c.f. Eq. (|24[) ): 
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ft)=/'^>)J4 (A9) 
x y \y J 

where 

Jq/M = J d PT fQ%(y,PT), (A10) 

and 

Qq(xq) = J d k T QQtt/q(xq, k T ). (All) 

The convolution formula ( |A9| ) can also be easily generalized to the case of n pions in the 
intermediate state, as in Fig.3(b): 

q n(S) = j^l... ^-J^j N {yx) q Q {y 2 /yi) ■ ■■qQ{yn/y n -i) qQ(x/y n ). (A12) 

J yi y-n 

As mentioned above, a formulation of the higher-order diagrams in terms of simple convo- 
lution integrals enables fewer integrations over coupled momenta to be performed, which 
simplifies considerably the numerical evaluation of these contributions. 



APPENDIX B: DIQUARK STRUCTURE 

In this Appendix we discuss DIS from possible diquark constituents in the nucleon. If 
one takes seriously the quark-diquark picture of the nucleon, in which the diquark itself is 
treated as an effective quasi-particle, then one must also consider the process where a diquark 
is struck by the DIS probe, Fig. 9, which then requires modeling the structure function of 
a bound diquark. (Although any realistic model of the nucleon, which is to reproduce the 
Callan-Gross relation, cannot have diquarks as elementary constituents.) For scalar diquarks 
this process was considered by Suzuki et al. ||14| . In a covariant framework scattering from 



S = 1 diquarks is more problematic, however, as was demonstrated in Ref. 0] for DIS off 
spin-1 meson configurations in the nucleon. Within the IMF approach discussed in Appendix 
A above, calculating the diquark distribution is considerably simpler — in fact one expects 
explicit symmetry between the constituent quark and diquark distribution functions: 
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wv) = /i%(i-i/). (bi) 



where fj)/^-(y) is the probability distribution of spin S diquarks in the nucleon with momen- 
tum fraction y. In a covariant treatment this symmetry, which is 'screened' by off-mass-shell 
effects, need not be explicit. 

The relation in Eq.(BT) is readily obtained from any QDN vertex function that is a 



function of the variable Mq D = (toq + Pt)/v + (w| + py)/(l — y) (see Eq.(^)) — for 



example a multipole form: ^ s > oc 1/(A| + Mq D ) us . With the interchange of rriQ *-> ms and 
y <-» 1 — y, this choice automatically results in the relation (|B1~1) |f23Tj3"4l , f40f . This constraint 
on the CQ-diquark symmetry consequently imposes constraints upon the shape of the quark 
distribution. Phenomenologically, the main difference between these forms arises at small 
x, where the Mq D - dependent functions give somewhat smaller distributions. The reason for 
this is the 1/y factor in Mq D , which at small x serves to suppress the quark distributions, 
which themselves depend on inverse powers of Mq D . 

Having outlined how one could systematically include DIS contributions from diquarks, 
we should point out however that DIS from diquarks is in fact a next-to-leading order 
process in the QQ coupling constant G. Hence it should be treated on the same footing as 
the order (l/f n ) 2 processes involving pion dressing, such as the direct scattering from the 
nucleon's pion cloud. While the higher-order (higher-mass) processes can in principle be 
systematically included, it is in fact more efficient to describe them in terms of the Regge 
model described in Section IV.B. 
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FIGURES 




FIG. 1. Leading twist contribution to deep- inelastic scattering of a photon (momentum q) from 
a constituent quark (p) within a nucleon (P), with \n) labeling the spectator quark system. 
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FIG. 2. DIS from a constituent quark with substructure, (a) direct photon diagram, which 
contributes to the quark distribution, (b) crossed photon diagram, associated with the antiquark 
distribution in the nucleon. 
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FIG. 3. Contribution to the nucleon structure function from the (a) low mass spectator (di- 
quark), (b) intermediate mass pion-exchange, and (c) large mass (asymptotic) components of the 
nucleon spectral function. 
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FIG. 4. (a) Total valence x{uy + dy) quark distribution, and (b) valence dy juy ratio, including 
only diquark spectators (process A), with two sets of mass parameters, mom = 1.0(1.2) GeV (solid) 



and m, 



0(1) 



0.6(0.8) GeV (dashed). The curves are evolved from i? = 1 GeV 2 to Q 2 = 10 GeV 2 , 



and the shaded area represents parametrizations of data taken from Ref. [ 30 1 . 
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FIG. 5. Typical dissociation process generated from the effective constituent quark-pion inter- 
action Cqtt. 
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FIG. 7. Ratio of the one-pion exchange (£>) to the handbag contributions (A) to the total 
valence distribution, q = uy + dy, for different Qtt vertex cut-off parameters. 
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FIG. 8. Total valence quark distribution evolved from n 2 = 1 GeV 2 to Q 2 = 10 GeV 2 , from 
process A and £>, and A alone (dashed), and from A — C together (solid). Shaded region represents 
data at Q 2 = 10 GeV 2 . 
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